NMR evidence for two-step phase-separation in Nd1.85Ceo.15Cu.O4_5. 
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By Cu NMR we studied the spin and charge structure in Nd2- x Ce x Cu04-<5. For x = 0.15, starting 
from a superconducting sample, the low temperature magnetic order in the sample reoxygenated 
under 1 bar oxygen at 900° C, reveals a peculiar modulation of the internal field, indicative for a 
phase characterized by large charge droplets ('Blob'-phase). By prolonged reoxygenation at 4 bar 
the blobs brake up and the spin structure changes to that of an ordered antiferromagnet (AF). 
We conclude that the superconductivity in the n-type systems competes with a genuine type I 
Mott-insulating state. 
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At present an important issue in cuprate superconduc- 
tivity is the nature of the electronic state(s) competing 
with the superconducting state (see Q, Q and references 
therein). In the hole doped cuprates evidence has been 
accumulating that at low dopings this competitor is a 
stripe phase [3j. It is believed that these stripe phases 
find their origin in the microscopic incompatibility be- 
tween the metallic- and Mott-insulating states jj, 0|. 
Stripes are a-priori not a unique way of resolving this 
incompatibility. This was sharply formulated recently in 
terms of the analogy with superconductivity stripes 
can be viewed as the analogue of the type II supercon- 
ductor while a type I behavior is also imaginable. In the 
latter the Mott-insulator and the metal/superconductor 
are immiscible and one expects a state characterized by 
droplets of metal/superconductor in a Mott-insulating 
background. The electron doped (n-type) superconduc- 
tors are microscopically quite different from the hole 
doped ones, while also macroscopic properties seem dis- 
tinct So far, no incommensurate spin fluctuations 
have been found by neutron scattering in this system 
[H E UH, arguing against stripe phases. Instead, re- 
cently it was reported that field induced commensurate 
antiferromagnetism re-emerges in the superconducting 
state in the presence of an Abrikosov vortex lattice [12j , 
suggesting that the superconductivity competes with a 
conventional antiferromagnet. 

Different from the p-type systems, the n-type systems 
stay insulating up to quite high dopings while the antifer- 
romagnetism of half filling degrades quite slowly. In fact, 
Vajk et al. |l0j showed that the behavior of this antifer- 
romagnet can be understood in detail assuming a random 
dilution of the Heisenberg quantum antiferromagnet, sug- 
gesting that the individual carriers stay strongly bound 
to impurity sites. Is the field induced antiferromagnet 
of a similar kind? Using NMR we will present here 
evidences that in between the superconductor and the 
strong pinning phase there is yet another phase. By re- 
oxygenating a superconducting Nd2- x Ce x Cu04 (NCCO) 



sample with x = 0.15 up to the point that supercon- 
ducting diamagnetism has disappeared we have detected 
a phase in a very small oxygen-doping range located in 
between the superconductor and the site diluted antifer- 
romagnet. This phase is extremely sensitive to further 
oxy genation. Approximately one oxygen atom per 100 
[13j unit cells suffices to destroy this phase allowing a 
site-dilution antiferromagnet to take over. With NMR 
one can measure the amplitude distribution of the mag- 
netic moments. These show a double-peak distribution 
in this special phase. Extensive Monte-Carlo simulations 
on a representative model show that such a distribution 
will only arise in either a stripe phase or in a situation 
where large droplets are formed. Since the neutron scat- 
tering appear to rule out the former we conclude that the 
superconductivity in the n-type systems competes with 
a genuine type I Mott-insulating state 0. 

Experimental. The experiments were performed on sin- 
gle crystals with x = 0.15, 0.12 and 0.08. For x = 0.15 the 
oxygen reduced sample (NCCOp) had a superconducting 
transition temperature of T c = 21 K. Half of the crys- 
tal (2x1x0.2 mm) was reoxygenated in air at 900° C till 
no superconductivity could be detected in the SQUID- 
susceptibility measurement (sample NCCOl). After per- 
forming the full set of NMR measurements, the same 
sample was oxygenated further under 4 bar of oxygen at 
850° C for 20 hours (sample NCC02), and NMR-scanned 
again. Sample preparation and characterization are de- 
scribed elsewhere [lj|. The Cu NMR spectra were ob- 
tained with conventional phase coherent pulsed NMR be- 
tween 4 and 350 K by sweeping the magnetic field at var- 
ious constant frequencies. The NMR line-shift was deter- 
mined by using a gyromagnetic ratio of 1.1285 MHz/kOe 
for 63 Cu and K — 0.238 % of metallic Cu as a reference. 

The Cu NMR spectra were studied for two orienta- 
tions (_? a ppi||c and B app \\\ab) of the crystals with respect 
to applied external magnetic field -B a ppi- An example for 
x = 0.15 with _? ap pi||c is shown in Fig. fjj, at high tern- 
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FIG. 1: Cu NMR spectra for x = 0.15 in NCCOp (for com- 
parison scaled down from 78 to 43 MHz) at 100 K (at this 
temperature the spectra for NCCOl and NCC02 are identi- 
cal to that of NCCOp) and NCCOl at 5 K. The broadening 
and shifts at low temperature are clearly visible. 



perature (100 K). The two narrow lines correspond to 
the central transitions of the two copper isotopes 63 Cu 
and 65 Cu and are superposed on a relatively broad back- 
ground which likely originates from quadrupole satellites 
with small quadrupole frequencies. For characterization 
of the samples, we measured the T dependencies of the 
normalized product of the line intensity (J) (corrected for 
T 2 effects) and T, the NMR line-sift 63 K, and the funda- 
mental relaxation probability W\ |l5|, [lT], [l8| . Here 
we concentrate on the low temperature field profiles on 
the Cu-sites only - the internal field effects are so large 
that the role of the Nd-moments can be neglected. 
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FIG. 2: (a) Amplitude distribution of the internal field for 
NCCOl at 5 K. The solid line gives the profile of the hyperfine 
field, as deduced from B \\ c. From the data for B J_ c the 
antiferromagnetic alignment of the electron spins in the (a, b)- 
plane appears to be strongly disordered, see text, (b) Simple 
reconstruction of the internal field distribution for a one and 
two dimensional spin structure. In ID the radial coordinate 
r with (incommensurate) wavevector k r have to be read as x 
and k x . In 2 D cylindrical symmetry is assumed. 

Analysis of the low temperature profiles. In the 
presence of the static hyperfine magnetic field -Bhf 



on the copper nuclei produced by ordered Cu 2+ 
spins, one should observe two Cu NMR lines at 
B = \ B l m \ + B li + 2S app iS hf cosa] a5 , where a is 
the angle between the directions of -B app i and Bhf. 
It is known from neutron scattering experiments on 
Ndi.85Ceo.i5Cu04 9] that the Cu 2+ spins are ordered 
antiferromagnetically in the afr-plane. Therefore, for 
B app i | |c (a = 90°) only one Cu NMR resonance is 
expected at a lower value than the field calculated from 
the gyromagnetic ratio, Bq. For those nuclei, where for 
some reason the internal field is cancelled, the resonance 
will be unshifted. The Cu NMR spectra recorded for 
B ap pi||c for x — 0.08, 0.13 and 0.15 can indeed be 
qualitatively explained in this way: for x = 0.15, see 
Fig. part of the nuclei mainly feel the external field 
(resonate at the fields of Fig. ^i) and the remainder 
experiences an internal field in addition. At lower doping 
the line shifts are larger. The spectra for -B app i||a6, see 
Fig. seem to be composed out of at least two lines 
of different intensities. This difference is not due to T 2 
or different 7r/2-conditions, because these are similar for 
both peaks (T 2 _1 = 23(1) x 10 3 s" 1 and tt/2 = 2 /is). 



For a quantitative analysis, the low-temperature cop- 
per NMR spectra of NCCOl at 5 K were simulated by 
exact diagonalization of the nuclear spin Hamiltonian Ti. 
with effective spin I = 1/2 (the quadrupolar interaction 
is weak compared to the magnetic interaction) for both 
isotopes 63 Cu and 65 Cu. For B \\ c, H can be written 
as Ti = 7^-B a ppi/ 2 + 7^-Bht [I x cos ip + I y sin tp] + ~fhB x I x , 
where the z-direction is along the crystallographic c-axis 
and ip is angle between -Bhf and the rf-field B\ (along 
the s-axis); for B J_ c, I z has to be replaced by I y . 
-B ap pi -L B\ in all experiments and the AF coupled elec- 
tron spins are assumed to be always in the (a, b) plane. 
The amplitude distribution of the internal field P(Bhf) 
was obtained via Monte-Carlo optimization to reproduce 
the data for B \\ c, see Fig. Analysis of the B _L c 
data (Fig. QJ;) with this pure in-plane AF spin model in- 
dicates that the alignment of the electron spins is rather 
disordered in the (a, 6)-plane, and on average show a field 
induced canting of some 20 degrees, which does not influ- 
ence significantly the distribution P(Bm) in Fig. In 
the x = 0.12,0.08 samples the computations give a sin- 
gle peak at a field ~ 3 T, consistent with a simple com- 
mensurate antiferromagnet. However, in the slightly oxy- 
genated x — 0.15 sample this spin-amplitude distribution 
is quite structured: the peak at 2.4 T (corresponding with 
a reduced moment of ~ 0.14/xs) is asymmetric showing a 
slow decrease on the smaller moment side, while a second 
peak at vanishing internal field corresponds with sites 
where the internal fields have completely disappeared. 

Interpretation. If the internal-field distribution would 
arise from a stripe-like spin structure with wave vec- 
tor k Xl we would have P(B) = d(k x x)/dB and hence 
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FIG. 3: 2D map of the blob (a) and stripe (b) results 
from thermalized averaged MC simulations, obtained with 
the hamiltonian discussed in the text. The spin-orientational 
disorder in the (a, 6)-plane is represented by a difference in 
color (complementary colors correspond to opposite aligned 
spins) , while the color intensity corresponds to the spin value. 
The spin-zero regions (black), as also seen in experiment, are 
found in the stripe simulations, but are almost absent for 
small blobs. The parameters for the blob and stripe panels 
are (J = 1, Vo = 0): Ji = 0.2, Q = 0.2, V s = 0.2, R = 10, 
T = 0.4 and J x = 1, Q = 0.4, V s = 0.2, R = 10 with T = 0.3. 




FIG. 4: Blob results from the MC simulations - a thermalized 
result for larger blobs than in the previous figure. Only if 
blobs are sufficiently large, zero field sites can be generated in 
the presence of the transferred hyperfine interaction. For still 
larger blobs the zero field contribution will further increase. 
Parameter values are J — 1, T — 0.3, Ji = 0, Q = 0.04, V s 
= 0.2 with R = 10 and Vc = - 1.5. 

x oc f Q B / Smax p( K B)dB. Using the experimental outcome 
for P(Bhf) the integration leads to the ID-wave depicted 
in Fig. ^jp; note that because the Cu-nuclei that experi- 
ence -Bhf have regular lattice positions, the wave vector 
needs to be incommensurate to realize such a field dis- 
tribution at the Cu-sites. The peak of P(Bm) at zero 
internal field in Fig. would then correspond with the 
nodes. However, the field profile for a radially symmetric 
two dimensional spin profile can be generated similarly 
via {k r rf cx / B/Smax P{B)dB and is given in the same 
figure. These simple calculations are only meant to show 



that both stripe like structures (comparable to the ID 
case) and charge blobs (2D case) are compatible with the 
experimental data. Furthermore the needed zero field re- 
gion in 2D is larger than in ID (due to the transferred 
hyperfine interaction the size of the blobs with zero spin 
exceeds the zero field region). 

The size enters into the problem even stronger 
because of the specific form of the hamiltonian. The 
internal field has not only a on-site contribution A, 
but also the transferred hyperfine interaction B w A/4, 
which gives a field of opposite sign and hence in a 
droplet a non-zero field will be generated on zero- 
spin sites (in a stripe configuration zero field states 
are easier to generate because it is an anti-domain 
boundary). To make this case more persuasive we 
performed Monte-Carlo simulations on a classical model 
designed to interpolate smoothly between stripe-like 
and droplet-like ground state textures [l^. This is a 
classical, spin-full lattice gas model which is similar to 
the model introduced by Stojkovic et al. [icf. It builds 
in the quantum-physical ingredient that electrons/holes 
('charge', ~ (1 — rii)) cause anti-phase boundaries in the 
spin system by a 'charge mediated' exchange interaction 
|2l| coupling the spins antiferromagnetically across 
the charge, H spin = J^2 x ^S x , y ■ (S x , y +i + S x+ i, y ) + 

^2x,y ^x.y • [(1 — n x , y+ i)S x .yJ r 2 + (1 — n x + i. y )S x +2,y] 

with J, Ji > 0. The frustration in the spin system 
is released by the formation of charged domain walls. 
The frustrated phase separation motive is build in by 
balancing the spin-mediated attractive charge-charge 
interactions with a repulsive 1/r Coulomb interaction 
of strength Q between the charges, which we cutoff 
at 10 sites (beyond the inter-stripe distance) for rea- 
sons of numerical efficiency 22]. We take a IF spin 
system and the spin-orientational disorder is build in 
by a spin pinning potential, breaking the ground-state 
degeneracy by choosing a direction (up to spin reflec- 
tion) for the antiferromagnetic background, as well as 
providing an inhomogeneous charge potential. It has 
the form H Pot = -V s \Y,x, y S *,y ■ P x ,y\, where P x , y = 

A^ye*®-'.*' eM-V(x-x') 2 + (y-y') 2 /R) 

with R the correlation length over which the antiferro- 
magnetic background is allowed to rotate. The angle 
Q x i, y i is chosen randomly with a uniform distribution 
and amplitude A X ' lV > — Random[0, 1]. Finally we also 
incorporate a quenched disorder potential acting on the 
charge in the form of Vc — Random [0, 1] 

Despite its simplicity this simple model is capable of 
generating quite complex ground states, ranging from 
highly disordered droplet patterns [J\ — > 0) to very or- 
derly stripe patterns J±, J >> Q,Vs. In the simula- 
tions we used small but finite temperatures and ther- 
mal annealing, mimicking the effects of quantum fluctu- 
ations in smearing the sharp textures |23j in this classical 
model. Two typical outcomes are illustrated in Figs. 13141 
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We have extensively scanned the parameter space of this 
model. Invariably, we find that a double peak structure 
is easily associated with stripe like patterns (Fig. 
but these configurations invariably lead to incommensu- 
rate scattering amplitudes, which are not observed exper- 
imentally. When droplets are formed instead, the zero- 
moment peak is always lacking for small droplets (Fig. [3] 
and Fig. 0J. Only droplets of 2 nm or more (25 sites 
or larger) mimic the peculiarities of the NMR spectra. 
This adds confidence to our claim that the double peak 
structure in the local moment distribution is the NMR 
fingerprint signalling the presence of large charge blobs. 
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FIG. 5: Copper NMR spectra (the intensity in arbitrary units 
along the vertical axis is plotted against the magnetic field 
value in Tesla) at 5 K for NCCOl and NCC02 (a). The lines 
in (b) are fits with an optimized internal field distribution, 
see text. 

How to reconcile our findings with the neutron scatter- 
ing evidences favoring the site-diluted antiferromagnet? 
The NMR sample has a volume of 1 mm 3 , much smaller 
than typical neutron targets and we found the large 
droplet phase by almost continuously monitoring the dis- 
appearance of the superconducting state while the sam- 
ple was oxygenated. Apparently, the large droplet phase 
in the NCCOl sample is linked to the presence of a tiny 
amount of apical oxygens (6 ~ 0.005 for NCCOl) 
The result for continuing the oxygenation of NCCOl in 
4 bar of oxygen as used in other studies is given in 
Fig. El The data show clearly that in NCC02 the phase 
of NCCOl is suppressed in favor of a new phase |24j . 
Applying a similar analysis as before, we find that this 
new phase is characteristic of an antiferromagnet with 
an internal field of 3 T - the phase seen in the neutron 
data. Although it is hard to arrive at precise numbers, 
it seems that one extra oxygen per roughly 100 unit cells 
suffices to fracture the droplets and to favor a state where 
the charge carriers are strongly bound to impurities. As 
the superconductor, the large droplet insulating phase 
shows an extreme sensitivity to oxygen concentration ex- 
plaining why this phase has been missed in the neutron 



scattering studies. 

In summary, we have presented evidence demonstrat- 
ing that the n-type cuprate superconductor is in a di- 
rect competition with a phase where the charge carriers 
form large droplets in a commensurate antiferromagnetic 
background. This phase is in turn extremely sensitive to 
the chemical conditions. Tiny amounts of excess oxygen 
suffice to destroy both superconductivity and the large 
droplet phase, in favor of a phase where the charge car- 
riers are strongly bound to the impurities. 
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